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Abstract 

Although the flux density map of a bulk superconductor provides in 
principle sufficient information for calculating the magnitude and the 
direction of the supercurrent flow, the inversion of the Biot-Savart law 
is ill conditioned for thick samples, thus rendering this method unsuit- 
able for state of the art bulk superconductors. If a thin (< 1 mm) slab 
is cut from the bulk, the inversion is reasonably well conditioned and 
the variation of the critical current density in the sample can be cal- 
culated with adequate spatial resolution. Therefore a novel procedure 
is employed, which exploits the symmetry of the problem and solves 
the equations non-iteratively, assuming a planar z-independent current 
density. The calculated current density at a certain position is found 
to depend on the magnetic induction. In this way the average field 
dependence of the critical current density J C (B) is obtained also at 
low fields, which is not accessible to magnetisation measurements due 
to the self-field of the sample. It is further shown that an evaluation 
of magnetisation loops, taking the self-field into account, results in a 
similar dependence in the field range accessible to this experiment. 
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1 Introduction 



At present, bulk superconductors with several centimeters in diameter and 
about one centimeter in thickness, trapping remanent magnetic fields ex- 
ceeding IT at 77 K, can be reproducibly grown pp. Texturing of the mono- 
lith is achieved by top seeded melt texture growth (TSMG), where crystalli- 
sation evolves from a seed crystal placed on top of the bulk. Five growth 
sectors are formed, propagating from the facets of the crystal through the 
entire material [2]. The critical currents achieved in each growth sector, 
especially as a function of the seed distance, are therefore of particular in- 
terest. 

A straightforward approach is to cut the sample and to characterise 
small pieces by magnetometry. However, this procedure is destructive and 
takes extensive measurement time. On the other hand, scanning Hall probe 
techniques can be employed to analyse the local properties of the bulk [3]. 
Among them the magnetoscan technique [S] proved to provide detailed in- 
formation on the critical current flow on a local scale. However, even the 
strongest permanent magnets used in the magnetoscan device activate cur- 
rents only in the uppermost layer of the bulk. The information refers to 
depths of less than 1 mm [B] and the inside of the bulk cannot be probed. 

Currents flowing in the entire volume of the bulk can be activated by 
performing a (zero) field cooled hysteresis loop in a magnet. In principle 
sufficient data to calculate the current density on the sub-mm length scale, 
where substantial changes in J c are expected, can be obtained from trapped 
field maps. 

Although elaborate and numerically stable techniques exist 00 [9], it 
was shown |10| that due to the large thickness of those samples the matrix 
inverted in such procedures is notoriously ill conditioned and can even be 
numerically singular. Thus, a reduction in thickness, either by grinding or 
preferably by cutting the bulk into disks, is mandatory to assess the local 
critical current distribution in this way. If the disk is sufficiently thin, a 
scanning mesh can be found which allows both an analysis of the current 
density on the sub-mm scale and a comparison to the magnetoscan signal 
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Figure 1: (Color Online) Estimate of the condition number for the param- 
eters used in the experiment (s = 0.2mm, d = 0.55mm). Since the matrix 
exceeds Scilab memory stack, the condition number was calculated for in- 
creasingly larger systems of up to 50 x 50 points. The convergence of the 
condition number is evident. (Inset) Symmetry of the matrix involved in the 
computation. The nine Toeplitz blocks constitute a block Toeplitz matrix. 

at certain positions. 

2 Numerics 
2.1 General 

In the coordinate system used in the following (cf. figure [4]) the top sample 
surface lies in the x, y-plane and the perpendicular component of the mag- 
netic induction B = B ■ e z is assessed. Similar to [7] the current density 
is expressed as the curl of a magnetisation density pointing in z-direction 
J = V x M e z . This implicitly satisfies current conservation V • J = for 
the z-independent planar current distribution J(x, y). Unlike in the basic 
equation of the problem is derived by splitting the magnetic induction into 
two components and using a scalar magnetic potential (see the Appendix 
for a detailed derivation) . Discretisation of the integral equation using cubic 
volume elements with constant M results in the 2D matrix equation 
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K i!jAl = F(A X ,Ay,Az) } »> » , , (2) 
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F(Ax,Ay,Az) = ^tan^f / AxAj/ =) . (3) 

1 y ' 4tt K Azy/Ax 2 + Ay 2 + Az 2 ' K ' 

Here c denotes the sample thickness, d the distance between the active 
area of the Hall probe and the top sample surface (gap), s the step width, 
and i, k = 1 . . . n x , j, I = 1 ■ • • n y the indices on the mesh; the antiderivative 
F is evaluated at the eight corners of the cubes in ([2]). 

Equation [1] can be mapped one-to-one to ID by substituting %' = i + 
j (n x — 1), k' = k + I (n y — 1). The equation now reads 

Tl x ' ' y 

^ K\ l>k ,M k , = B V (4) 

k'=l 

and the problem can be tackled by matrix inversion algorithms. 

As pointed out earlier, the matrix Kij^i is in fact a Toeplitz block 
Toeplitz matrix [7] . This is a consequence of the translation invariance of the 
Biot-Savart law (Ki t j^,i = ^|i-fc|,o,b'-Z|,o) an d results in a highly symmetric 
matrix K' if (cf. inset of figure [T]). Therefore, both efficient storage and a 
fast algorithm for solving the system can be expected. However, [7] exploits 
the symmetry only for the storage of the matrix elements and the method of 
conjugated gradients using the fast Fourier transform (FFT) is employed to 
solve the linear equation. Although the procedure is fast, employing FFT 
[7J [8] implies the unnecessary assumption of periodicity in B outside the 
measurement area, which may create artefacts, if currents are flowing close 
to the edge of the scanning area. 

Block Toeplitz matrices occur in a number of problems, such as image 
reconstruction or system identification. Fortunately, an efficient and fast 
algorithm, which exploits the symmetry of the structured matrix, is provided 
in [11]. The computation time is approximately 10 minutes for a 200 x 
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200 system on a desktop PC. This is presumably longer than FFT based 
algorithms, but still much less than the actual measurement time. 



2.2 Condition number 

The residuum of the solution is of the order of the machine error. However, 
it is shown in standard numerical algebra textbooks, that matrix inversions 
can amplify a relative (measurement) error in the right-hand side of @, 
leading to an unknown error in the calculated magnetisation density cm- 
This behaviour is described by the condition number k of a linear system 



where the errors are measured in Euclidean norm. It was shown in [10] 
that the condition number of the inverted matrix is in notoriously high for 
bulk superconductors. Therefore, special attention has to be paid to the 
choice of the step width s, defining n(K') at a given gap d. With a real- 
istic distance d of about 0.15 mm (see section [3|) and a sample thickness of 
0.55 mm, the condition number was estimated using the numerical compu- 
tation package Scilab [12] (cf. Fig. [[]). 

As a rule of thumb one aims at a condition number k such that 



holds. It was found that the matrix inversion is reasonably well condi- 
tioned as long as the step width is larger than the gap. The actual choice 
of s = 0.2 mm is adequate to resolve the spatial variations of the critical 
current on a sub-mm scale. The condition number k ~ 9 combined with 
a rough estimate of the relative measurement error 1 1 1 1 2 ~ 0.01 results in 
IkAflb ~ 0.09, satisfying (jSJ). 

Calculating the current density involves computing a (numerical) deriva- 
tive, which implies that the relative error in J(x, y) will be high, if the change 
in M(x, y) (the current at this position) vanishes, as for example outside of 
the bulk or in defects. This is a peculiarity of measuring the relative error of 
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a quantity close to zero. Note that the absolute error in J{x,y) is bounded 
by © and expected to be acceptably low for most of the current distribution 
inside the sample as long as ([6]) holds. 

3 Experimental 

A thin disk was cut from an undoped YBCO bulk superconductor with 
a diameter of 26.5 mm, which was grown by the top seeded melt growth 
technique [3], using a diamond saw. The cut was made near the upper 
surface of the bulk and the disk was polished with abrasive paper to a 
thickness of 0.55 mm. 

Commercial Hall probes from Arepoc with active areas of 25 x 25 /im 2 
(trapped field map) and 50 x 50 /im 2 (magnetoscan) were used in the ex- 
periments. All scans were carried out with a step width of 0.2 mm, which 
allows to apply the inversion (see section I2.2p . 

For the trapped field maps a scanning area of 28 x 28 mm was used. 
The sample was field cooled in a split coil magnet in a field of 1.4 T. In 
order to minimise relaxation effects during the measurement, the scan was 
started 10 min after sweeping the magnetic field to zero. The temperature 
of the liquid nitrogen bath was recorded prior to and after the scans and was 
found to be 77.2 K, increasing due to oxygen uptake by about 0.1K during 
the measurement. The resistive offset of the Hall probe was determined at 
the end of each measurement, with the Hall probe still immersed in liquid 
nitrogen, but at a large distance from the bulk. 

A somewhat wider scanning area of 30 x 30 mm was used for the mag- 
netoscans. This precaution was taken to assure that the magnet would not 
produce artefacts by stopping over the bulk after finishing a single line of 
the scan. The magnetoscan was carried out using a small SmCo permanent 
magnet of 2 mm in diameter, applying an induction of around 150 mT at 
the top sample surface of the bulk. Further details of the technique can be 
found in [5]. 
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4 Results 



4.1 Inversion of the Trapped Flux Density Profile 

Assuming that the variation of the critical current over the sample thickness 
can be neglected immediately implies that the trapped field profiles recorded 
on the bottom and top surface are identical, which was confirmed by the 
experiment. The maxima of the trapped field were equal within one percent 
on both sides of the bulk and found to be 252 mT. Strong negative fields of 
up to — 100 mT were detected close to the edge of the sample, which result 
from the high diameter to thickness aspect ratio of the disk (cf. Fig. [2^). 

The high reproducibility between several measurements shows that the 
gap between the Hall probe and the sample surface (~ 0.15 mm) remains 
unchanged in subsequent runs. 

The inversion of the trapped flux density profile of the bottom surface of 
the disk is depicted in figure [5b. Most of the bulk carries a current density 
of around 4 • 10 8 Am -2 in the remanent state. The defects (blue arrows) 
close to the edge appear as regions, where the critical current is drastically 
reduced, most likely due to cracks or large scale inclusions. A remarkable 
result of the inversion is the detection of the c-growth sector at the center of 
the bulk. It is clearly displayed in the current density map as a rectangular 
area with low critical current density (green region). 

In addition, a clear negative correlation between the critical current den- 
sity and the magnitude of the perpendicular magnetic induction B is found, 
for example close to one of the a— a growth sector boundary, where both the 
magnetic field and the current density change simultaneously (small white 
arrows). The correlation is most prominent at low fields, i.e. high current 
densities of up to 10 9 Am -2 flow close to the sample edge, where the mag- 
netic induction changes sign and therefore B = 0. Moreover, a point inside 
the bulk (lowest white arrow) with reduced magnetic induction is repro- 
ducibly detected, where the current density significantly exceeds the nearby 
current densities. This demonstrates that a strong field dependence of the 
critical current is present, especially at low fields. 
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(a) Flux density map 




Figure 2: (Color online) (a) A number of extended defects are detected in 
the trapped flux density map (blue arrows) and strong negative fields of up 
to — 100 mT are observed (green area), (b) The inversion reveals a strong 
correlation between the magnetic induction and the current density (white 
arrows in a, b). In the defects (blue arrows) as well as in the c-axis growth 
sector (green rectangular area in the center) the critical current density is 
significantly reduced, (c) The c-axis growth sector is again detected by 
the magnetoscan. The small length scale of the variations in the signal 
indicates strong local changes in the critical current (gray arrow). (The 
black rectangles indicate the cubes cut for magnetometry.) 
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4.2 Comparison to Magnetoscan 

Due to the strong field dependence of the critical current, it is difficult to 
compare different regions of the bulk, as the self-field of the bulk is position- 
dependent in the remanent state. Contrary, in the magnetoscan the back- 
ground field of the permanent magnet is constant and the self-field is smaller, 
since the currents are activated only in an area of about the magnet's diam- 
eter [6]. 

All except one of the prominent defects are found in the magnetoscan 
(cf. Fig. [2b) . A possible explanation would be a large defect or inclusion at a 
depth, which exceeds the penetration depth of the permanent magnet. Sim- 
ilar to the inversion, the reduction of the critical current is evident from the 
low magnetic response in the c-axis growth sector. Moreover, the observable 
granularity on a sub-mm scale of the bulk's response indicates strong local 
variations in the critical current density. 

4.3 Comparison to Magnet ometry 

The calculated current densities can be correlated with the z-component 
of the magnetic induction at their position. For this purpose the magnetic 
induction in the central plane inside the bulk was calculated after the in- 
version (see the Appendix), resulting in a J C {B) plot (cf. figure [3]). Only 
points with positive induction B > were considered, which effectively cuts 
off the noise due to the extended defects close to the sample edge. 

For comparison four cubes (2 x 2 x 0.55 mm 3 ) were cut from different 
positions of the bulk (cf. Fig. [2b) and analysed in a SQUID magnetometer 
at the temperature of the liquid nitrogen bath during the previous scans. 
The magnetisation loops were evaluated assuming the Bean model. A large 
sample to sample variation was found, at this point the lowest critical cur- 
rent density being located at the c-axis growth sector in the center of the 
bulk (cubes #3 and #4), in agreement with the above results. Although 
numerical errors in the inversion procedure (mainly due to the numerical 
derivation) cannot completely be excluded, the variation in the magnetisa- 
tion data elucidates the scatter in the calculated current densities, if one 
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Figure 3: (a) J C {B > 0) correlation from the inversion, showing the average 
field dependence of the critical current density. The magnetisation J C {H) 
data is similar to the smoothed data from the current calculation (green 
line) for high fields. The curves differ at low fields, because the self-field 
is not taken into account in the J C (H) data of the cubes. Further, the 
critical currents vary also from sample to sample in the magnetisation data; 
the cubes cut from the seed area show the lowest critical current density. 
(6) Evaluation of the hysteresis loop for cube #1 taking the self-field into 
account. The field generated by the supercurrents in the sample shifts the 
applied magnetic induction to higher fields (arrow), resulting in a J C {B) 
dependence, which is now similar to the average data obtained from the 
current calculation. There is no indication for a flattening of the curve at 
small fields. 
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takes into account that the averaged volume is much larger in the SQUID 
measurements. This suggests that the scatter is primarily due to strong lo- 
cal variations of the critical current in the bulk, which is further supported 
by the granularity observed in the magnetoscan signal (cf . Fig. [2b) . 

The average field dependence of the critical current density from the 
magnetisation loops J C {H) shows good agreement with J C (B) obtained from 
the trapped flux density maps, except at low fields, where the curves clearly 
differ (cf. Fig. (3^i). This can be explained in terms of the sample's self-field, 
which is not considered in the evaluation of the magnetisation experiment, 
since it takes only the externally applied field /j,qH into account and neglects 
the field generated by the currents flowing in the sample. 

To estimate this contribution the magnetoscan signal serves as a first 
approximation, because the currents are induced in an area determined by 
the magnet's diameter and the volume of current flow mainly contributing 
to the magnetoscan signal is therefore similar to the cubes used in mag- 
netometry. Indeed the average induction of the magnetoscan (~ 70 mT) is 
close to the field, where the J C {H) and J C (B) start to differ. 

For a more detailed analysis an evaluation method was applied (again 
based on the Bean model), which accounts for both the externally applied 
field and the mean self-field in the sample, thus providing an approximate 
J C (B) dependence [13]. Especially for low applied fields, where self-field 
effects become dominant, the evaluation shifts all data points to higher 
fields (cf. Fig. [3b). This effect is particularly clear for the remanent state 
at zero applied field (arrow), where the magnetic induction is solely due 
to the trapped self-field. Consequently, a sample cannot be probed at zero 
magnetic induction by SQUID magnetisation experiments. 

When accounting for the self-field contribution, good agreement between 
the J C (B) curve from the SQUID loops and the average J C (B) curve from the 
inversion is found in the range accessible to the magnetisation experiment 
(depicted in figure [3b for one of the cubes). Moreover, the field, where the 
correction starts to become effective and the self-field becomes important, 
is equal to the simple estimate made above using the magnetoscan data 
(70 mT). 
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In contrast to SQUID magnetometry, the current calculation allows one 
to analyse the J C (B) dependence in fields ranging from zero induction to the 
maximum trapped field. A clear dependence of the critical current density 
on the magnetic induction B is revealed in this way. There is no indication 
for a flattening in the average J C {B) at low fields. 

5 Summary 

The inversion of the Biot-Savart law represents an ill conditioned problem 
for bulk superconductors, but parameters, e.g. the step width of scan and 
the thickness of the sample, can be found, which allow its application to 
thin disks cut from the sample. The matrix equation was solved without 
any additional assumptions by a fast algorithm, which exploits the symme- 
try of the problem. In this way, for example the c-axis growth sector was 
clearly identified as a distinct region of low critical current density, a result, 
which is also obtained by the magnetoscan and confirmed by magnetisation 
measurements. 

Prom a technological point of view both methods, the magnetoscan and 
the inversion of the flux density map, are complementary. The magnetoscan 
provides information on the local critical current density at a certain con- 
stant background field, which is of interest for superconducting bearings. On 
the other hand, information about the remanent current flow at the bulk's 
self-field is important for magnet applications and can be obtained by the 
inversion of the trapped flux density map of thin disks cut from the sample. 

In addition, the current calculation was found to provide the important 
average field dependence of the critical current also at low fields (below 
the self-field), a region, which is not accessible to magnetic measurements, 
even if the self-field is explicitly accounted for in the evaluation procedure. 
There is no indication of a plateau in J C {B) and the critical current density is 
found to depend on the field in a continuous way also at the lowest magnetic 
inductions. 
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Figure 4: Coordinate system employed. 
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Appendix 

The magnetostatic Maxwell equations (V • B = and Ampere's law Vx5 = 
Ho J ) are treated by splitting the magnetic induction into a sum of two fields 

B = Vo(n + M), (A-l) 
where M is chosen to satisfy the differential equation 



V x M = J (A-2) 

in entire space, which satisfies current conservation V • J = 0. Equa- 
tion (|A-2p defines M apart from a gradient field and an arbitrary constant, 
which are both chosen to be zero. Consequently M vanishes outside the 
sample and can therefore be interpreted as a magnetisation, which is con- 
fined to the sample volume and establishes the current density J by spatial 
variations. Ampere's law results now in a homogeneous equation for the 
field fl 
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Vxfi = J-VxM = 0, 
which can thus be derived from a scalar potential 



(A-3) 



- V$ = O . (A-4) 

Taking the divergence of this expression and substituting the second 
Maxwell equation (V • B = V • (O + M) = 0) results in 

V 2 $ = V • M (A-5) 
and the problem can be solved by the method of Green's functions 

,p) = -2. /V^^O (A . 6) 

47r 7 |r-r'| 47r / \r-r'\ 

The planar z-independent current flow can be represented by M = 
M(x,y) e z and consequently the first term vanishes since V • M = d z M z = 
in the sample volume. Further, only the surfaces perpendicular to M con- 
tribute to the second term. Using (|A-6|) the measured induction B is ex- 
pressed solely function of M 



^^^y/^^W (A - 7) 



</ 



Q(r) 



Here Ax = x — x' , Ay = y — y', Az = z — z' and the notation for 
evaluating antiderivatives are used to indicate the positive contribution from 
the top (V = o?) and the negative from the bottom {z' = d + c) surface. 

Note, that ft vanishes for the case of an infinitely long slab (zero de- 
magnetisation) as Az — > oo in (IXTj) . In this case the induction is simply 
determined by the magnitude of M at a certain position and B = [IqM 
inside and B = outside the slab, where M = 0. However, for finite geome- 
tries there will always be a contribution from Q in (|A-7p and the relation 
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between B and M is non-local. Therefore the integral equation (|A-7p must 
be solved in order to attain M. 

The corresponding set of linear equations for the discrete measured data 
is formulated by approximating the sample as an array of cubes with con- 
stant M. Summing over all elements results in the matrix equation 

n x n y 

K i,j,k,iM kt i = /i(Aj = Bij . (A-8) 

k=l 1=1 

Here, the matrix entries Ki j k i are calculated by evaluating the an- 
tiderivative 



F(Ax,Ay,Az) = £ 
Mo 



dx dy 



Az 



— tan 1 (- 



(Ax 2 + Ay 2 + Az 2 ) 3 / 2 
AxAy 



4vr " Azy/Ax 2 + Ay 2 + Az 2 ' 
at the eight corners of the cubes 



(A-9) 
(A-10) 



K 



i,j,k,l 



F(Ax,Ay,Az) 



s(i-fc+|) 
s(i-fc-i) 



d+c 



(A-ll) 



Once the components M^i are obtained by inverting (|A-8|) . the current 
density can be calculated by employing (|A-2p . Further, the induction at any 
distance d outside the sample can be obtained from (|A-8p . Note, that the 
magnetic induction B in the central plane (d = — c/2) of the sample, which 
is needed to obtain J C (B), is given by 



n x fly 

HoMij + ^2^2Kij t k,iMk,i = fio(Mij + ftij) = Bij 

k=l 1=1 

as M does not vanish inside the sample. 



(A-12) 
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